Charge density wave fluctuations, heavy electrons, and superconductivity in KM2S2 
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Understanding the complexities of electronic and magnetic ground states in solids is one of the main goals 
of solid-state physics. Materials with the canonical ThCr2Si2-type structure have proved particularly fruitful 
in this regards, as they exhibit a wide range of technologically advantageous physical properties described by 
"many-body physics," including high-temperature superconductivity and heavy fermion behavior. Here, using 
high-resolution synchrotron X-ray diffraction and time-of-flight neutron scattering, we show that the isostruc- 
tural mixed valence compound, KN12S2, displays a number of highly unusual structural transitions, most no- 
tably the presence of charge density wave fluctuations that disappear on cooling. This behavior occurs without 
magnetic or charge order, in contrast to expectations based on all other known materials. Furthermore, the 
low-temperature electronic state of KM2S2 is found to exhibit many characteristics of heavy-fermion behavior, 
including a heavy electron state (m* /m e ~ 24), with a negative coefficient of thermal expansion, and supercon- 
ductivity below T c = 0.46(2) K. In the potassium nickel sulfide, these behaviors arise in the absence of localized 
magnetism, and instead appear to originate in proximity to charge order. 



I. INTRODUCTION 

Quantum coherence of electronic states in metals, or more 
generally "many-body," emergent phenomena, such as su- 
perconductivity, result from electron-electron or electron- 
phonon interactions established by the constraint of the lat- 
tice. Many materials that give rise to many-body physics 
(i.e., high-temperature superconductivity or heavy-fermion 
behavior) are comprised of layers of edge- sharing [MX4] 
tetrahedra, where M is a transition metal and X is a main- 
group element, as commonly found in the Th&2Si2 (e.g., 
Bai_ x K x Fe 2 As2, ffl K x Fe 2 _ J Se 2 P or URu 2 Si 2 E] ) or ZrCuSiAs 
(e.g., SmFeAsOi_ x Fj31) structure-types. Extensive work has 
found rich electronic phenomena in these materials, includ- 
ing hidden-order in URu2Si2, 5 nematic orderP^J valley density 
wave orderpl magnetoelastic coupling, 9 11 and Fermi surface 
nesting. 12 It is a generally accepted fact that the presence of 
magnetism and/or magnetic fluctuations are important in pro- 
ducing the correlated electronic behavior in these materials. 

Here, we report that KM2S2 has a similarly rich struc- 
tural and electronic phase diagram in the absence of local- 
ized magnetism, with several features unexpected under tra- 
ditional theories of strong electron interactions including (1) 
the disappearance of charge density wave (CDW) fluctuations 
concomitant with an unusual increase in local symmetry on 
cooling without trivial charge order, and (2) an enhancement 
of the effective conduction electron mass at low temperatures 
coupled with negative thermal expansion. The former is un- 
expected on thermodynamic grounds as the increase in local 
symmetry implies a decrease in the configurational entropy 



of the structure. The latter, an increase in electronic entropy, 
is a hallmark of the many-body "heavy-fermion" state, but is 
unexpected as KM2S2 shows no signs of the localized mag- 
netism associated with producing such a stated 

Instead, our findings are most consistent with KN12S2 har- 
boring electronically driven phase transitions that arise from 
changes in hybridization of a bath of delocalized conduc- 
tion electrons with localized and bonded (i.e., CDW) elec- 
trons making it an ideal compound for study of the coupling 
between charge and structural degrees of freedom in mixed- 
valence materials. Furthermore, these results demonstrate that 
proximity to charge order alone, without localized magnetism, 
can drive strongly correlated physics and warrants further ex- 
perimental and theoretical attention. 



II. METHODS 

Polycrystalline, lustrous and orange-yellow powder of 
KN12S2 was prepared as previously described, but with a sub- 
stitution of S for Se.^J All samples were prepared and han- 
dled exclusively inside an argon-filled glovebox; no impu- 
rities were detected by laboratory X-ray diffraction. High- 
resolution synchrotron X-ray diffraction data were collected 
using the high-resolution powder diffractometer at the Ad- 
vanced Photon Source on beamline 11 -BM 15 from polycrys- 
talline powders sealed in an evacuated fused silica capil- 
lary backfilled with /? He = 10 torr. Data for T < 100 K 
were collected using a He cryostat (Oxford Instruments); 
data collected at T > 100 K were collected using a ni- 
trogen cryostream. Rietveld analyses were performed us- 
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The synchrotron X-ray diffraction 
(SXRD) data revealed the presence of 3 wt% K2M3S4 impu- 
rity and a 1 wt% N13S2 impurity, which were included in the 
Rietveld analyses. The chemical occupancies were refined in 
Rietveld analysis of the SXRD data to test if the KM2S2 phase 
was substoichiometric; the values varied less than 1% from 1, 
thus the values were fixed to unity. 

Neutron total scattering data were collected at temperatures 
between 5 K and 300 K on poly crystalline KN12S2 (loaded in 
a vanadium can with a He atmosphere) using the time-of flight 
HIPD and NPDF instruments at the Lujan Center, Los Alamos 
Neutron Science Center, Los Alamos National Laboratory. 
PDF analysis was performed on the total neutron scattering 
data and G(r) were extracted with g max = 29 A -1 (HIPD) and 
Gmax = 35 A -1 (NPDF) using PDFgetNP Least-squares fits 
to the PDF were performed using PDFgui. 19 Reverse Monte 
Carlo simulations of 20x20x6 supercells (24000 atoms, ~75 
A/side) were performed using RMCprofile, 20 while applying 
a small penalty for breaking tetrahedral coordination and a 
hard- sphere cut off for the Ni-S bond distance. The structural 
parameters (Ni-Ni distance and bond valence sums) were 
compiled from the average of four independent simulations. 
Atomistic visualization was accomplished using VESTAL 

Physical properties were measured using a Physical Prop- 
erties Measurement System, Quantum Design, Inc; for mea- 
surement below 1.8 K, a dilution refrigerator option was used. 
Specific heat measurements were performed using the quasi- 
adiabatic heat-pulse technique on sintered poly crystalline pel- 
lets attached to the sample stage using thermal grease. Mag- 
netization measurements were carried out at juoH — 1 T and 
2 T, and with the susceptibility estimated as % « AM / AH = 
[Mxr — Mit]/[1T]. Isothermal, field-dependent magnetization 
measurements were performed over a range of temperatures 
and at fields from juoH = to 9 T to determine the fraction 
of impurity spins that contribute to % at low-temperatures. 
A self-consistent global fit of the field-dependent magnetiza- 
tion data sets was performed to a Brillouin function for impu- 
rity paramagnetic spins (with a single set of three parameters 
g = 2, /, and concentration, fixed to be the same at all temper- 
atures) and to a linear function (the true temperature depen- 
dent susceptibility). For resistivity measurements, platinum 
wires were attached to sintered polycrystalline pellets using 
silver paste and dried under argon in a four point configura- 
tion. Equivalent results were achieved by using Gao.85Ino.15 
as a molten solder. 



III. RESULTS 

A. Synchrotron X-ray and Time-of-Flight Neutron Scattering 

Analysis of high-resolution synchrotron X-ray diffraction 
(SXRD) data indicates that the average crystallographic sym- 
metry of KM2S2 is tetragonal (IA/mmm) at all temperatures 
measured (6.4 K < T < 440 K; Figure [TJ. In KNi 2 S 2 , the 
[Ni2S2]~ layers of edge-sharing [MS4] tetrahedra are sep- 
arated by K + ions [Figure [2) ; this leaves the nickel atoms 
with a formal valence of "Ni L5+ ". No periodic distortions 
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FIG. 1 : Representative Rietveld analyses of SXRD data collected at 
(a) T = 300 K and (b) T = 6.4 K . 



TABLE I: Structure parameters of KM2S2 obtained from Rietveld 
refinement of data collected from synchrotron X-ray powder diffrac- 
tion, showing fractional coordinates (x,y,z), chemical site occupan- 
cies (occ), and isotropic atomic displacement parameters (Ui SO ) de- 
scribed with the IA/mmm (139) spacegroup. 
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are found, in contrast to the commensurate distortions ob- 
served in KCu2Se2 22 or the incommensurate modulation of 
P-SrRh 2 As 2 P 

Nonetheless, detailed analysis of the diffraction data reveals 
a structural transition near T ~ 75 K and negative thermal ex- 
pansion below r~9K. The temperature dependence of the 
unit cell volume, extracted from Rietveld analysis, is shown in 
Figure [2] A 1 wt% impurity (N13S2, undetectable by labora- 
tory XRD) included in the SXRD data refinements acts as an 
internal standard. First, there is clear change in slope near T 
= 75 K. Second, the unit cell volume remains constant from T 
= 13 to 8.9 K, but then increases with further cooling [Figure 
|2jb)]. The former is consistent with a structural change at T = 
75 K, 24 while the latter indicates a switch to negative thermal 
expansion behavior. 

The change in slope in unit cell volume at T = 75 K is cou- 
pled to a striking change in the average structural model re- 
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FIG. 2: (a) From analysis of temperature-dependent high-resolution 
synchrotron X-ray diffraction (11 -BM, Advanced Photon Source), a 
clear transition in unit cell volume of KNI2S2 is observed at T ~ 
75 K (dot-dashed line); no transition is observed at T = 250 K (dot- 
ted line), (b) The temperature dependence of the unit cell volume 
of KN12S2 shows two anomalies below T < 13 K, with an overall 
negative coefficient of thermal expansion below I < 8.9 K. (dashed 
lines guide the eye; inset: unit cell of KM2S2). 



quired to describe the data. Adequate fits to the SXRD data 
collected at room temperature are obtained only when the Ni 
atoms are displaced off of the high symmetry position and are 
instead statistically distributed in a split- site model with the 
lower 2mm site symmetry on the 8g (0, 0.5, z) Wycoff position 
[Figure [TJa), Table [I). The strongest evidence for the presence 
of this distortion comes directly from a Fourier difference map 
generated from the undistorted structural model, illustrated in 
Figure [3] A dearth of scattering intensity is located above and 
below the 4d Ni site, while an excess of intensity is located 
at the Ni position of (0, 0.5, 0.25), illustrated in Figure |3jb). 
Introduction of the distortion produces an undisturbed Fourier 
difference map [Figure [3jc)] . In contrast, at T = 6.4 K, the 
SXRD data are described by an ideal ThCr2Si2 structure, with 
Ni atoms on the Ad (0, 0.5, 0.25) Wycoff position with 4m2 
site symmetry [Figure [TJb), Table [I). The temperature depen- 
dence of the off-centering is shown in Figure |3jd). There is 
only a small variation with temperature above T > 75 K. How- 
ever, the off-centering abruptly disappears on cooling below 
T = 75 K. 

Since there are no supercell reflections corresponding to a 
long range periodic order of the off-centering, pair distribu- 
tion function (PDF) analysis of neutron total scattering data 
was used to probe the nature and spatial extent of the dis- 
tortions. Figure |4ja) shows the PDF analysis of total scat- 
tering data collected at T = 300 K. Consistent with the off- 
centering, there are significant shoulders to the peak at r ~ 
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FIG. 3: (a) Illustration of the (100) crystal plane of KNi 2 Se 2 and 
Fourier difference maps from data collected at T = 300 K, viewed 
on the (001) plane reveal the displacement of Ni atoms from the (b) 
ideal 4d (0, 0.5, 0.25) position to the (c) displaced site at 8g (0, 0.5, z). 
The schematic unit cell on the right illustrates the location of the Ni 
atoms (blue spheres), (d) Coincident with the transition in unit cell 
volume are charge density wave (CDW) fluctuations that disappears 
on cooling as T < 75 K; these CDW fluctuations manifest as a dis- 
placement of the Ni atom (z) from the ideal 4d (0, 0.5, 0.25) site to 
the 8g (0, 0.5, z) site. Through all measured temperatures, the space- 
group symmetry, lA/mmm is retained, with the displaced z position 
requiring a split- site occupancy (4 d to 8g). The color denotes the 
normalized Rietveld goodness-of-fit (% 2 /%best ; brighter = better fit) 
to show the significance of the displacement shift at T ~ 75 K. 



2.68 A corresponding to modulations in nearest-neighbor Ni- 
Ni distances. These displacements again only have a weak 
temperature-dependence at high temperature, but abruptly dis- 
appear below T ~ 75 K, concomitant with the observations 
from the SXRD analysis. Further, the ideal lA/mmm crys- 
tal structure from the SXRD analysis provides an excellent fit 
to the T = 15 K PDF [Figure [Jc), LS]. Ripples in the PDF 
with a period Ar = 2n/Q mSiX = 0.18 A, amplified at low r, 
are artifacts the finite Fourier transformation used extract the 
PDF from the scattering data. 25 In contrast, while a split-site 
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FIG. 4: (a) Temperature dependence of the neutron pair distribu- 
tion function derived from neutron total scattering data (HIPD, Lujan 
Center) of KNI2S2 shows the presence of multiple nearest-neighbor 
Ni-Ni distances observed at T = 300 K that disappear on cooling 
(dashed line: ideal Ni-Ni separation), (b) Quantitative modeling 
of the PDF (NPDF, Lujan Center) measured at T = 300 K reveals 
a poor least-squares fit of the IA/mmm crystal structure determined 
by SXRD (LS) due to the presence of multiple Ni-Ni distances (ar- 
rows). Reverse Monte Carlo (RMC) simulations accurately describe 
the data, (c) In contrast, the T = 15 K PDF (NPDF) is quantitatively 
described by a least-squares fit of the average crystallographic struc- 
ture with a single Ni-Ni distance. 



displacement of the Ni position describes the SXRD data at 
T = 300 K, it does not adequately describe the shoulders of 
the nearest-neighbor Ni-Ni correlation peak in the T = 300 K 
PDF [Figure |4jb), LS and arrows]. 

The T = 300 K PDF requires at least three distinct Ni- 
Ni distances at r ~ 2.57(1), 2.70(1), and 2.86(1) A (Fig- 
ure [5]), which are not provided by the split- site model used 
for Rietveld analysis of the SXRD data. Deconvolution of the 
nearest-neighbor pair- wise correlations is provided by fitting a 
linear combination of tetragonal (IA/mmm) and orthorhombic 
(Fmmm) phases to the PDF (Figure [5]). Each phase has split- 
site occupancy of the Ni atoms which are displaced along the c 
axis. All non- special internal coordinates and unit cell dimen- 
sions were allowed to refine along with the relative contribu- 
tion of each phase. The resulting relative contribution of each 
fraction of each phase is fu/ mmm = 48 at% and fFmmm = 52 
at% from least- squares refinement. From the models, we ex- 
tract three nominal Ni-Ni distances: ^Ni-Ni = 2.56, 2.67, and 
2.85 A, as illustrated at the right of Figure [5] However, this 
structural model only describes the first coordination sphere 
of the Ni-S and Ni-Ni correlations; beyond r > 3.5 A, even 
this multi-phase model fails to describe the PDF. 

Reverse Monte Carlo (RMC) simulations of the neutron 
pair distribution function and Bragg profile produce atomistic 
configurations that are compatible with the average crystal- 
lographic symmetry and extended pairwise correlations. [Fig- 
ureQb), RMC and Figurega)] . Projection of all 24000 atoms 
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FIG. 5: Least-squares fit of a linear combination of tetragonal 
IA/mmm and orthorhombic Fmmm models of KN12S2 to the T= 
300K PDF allows extraction of the local Ni-S and Ni-Ni distances; 
however, the structural model fails to describe the data beyond r > 
3.5 A. 



from the large supercell back onto the crystallographic unit 
cell resembles the anisotropic atomic displacement parame- 
ters obtained from Rietveld analysis [Figure [6jb)]. Further- 
more, the supercell does not reveal any locally ordered pat- 
terns [Figure [6|c)] . Instead, short and long Ni-Ni bonds ap- 
pear randomly distributed throughout the lattice [Figure [6|d)]. 

Statistical analysis of the resulting RMC supercell yields 
an equivalent ensemble distribution of Ni-Ni distances, while 
also describing the extended pair- wise correlations. The tri- 
modal histogram of Ni-Ni distances [Figure |7ja)] appears 
with maxima centered around r ~ 2.56, 2.67, and 2.85 A, 
consistent with the least-squares analysis [Figure [5]. To en- 
sure this distribution is not the trivial result of harmonic but 
anisotropic atomic displacements, artificial PDFs and Bragg 
profiles of KM2S2 were generated using an ideal IA/mmm 
unit cell with anisotropic thermal displacements, akin to 
Ref . [26] These profiles were fit using RMC simulations with 
the same starting supercell as used with the experimental data. 
The histogram of Ni-Ni displacements from an anisotropic, 
but harmonically distorted structure has much more symmet- 
ric and singly distributed peak shape. Subtraction of the con- 
trol simulation from experimental distribution emphasizes and 
confirms the presence of three populations of bond lengths, as 
illustrated in Figure |5ja). From the RMC supercell, we were 
able to extract an ensemble of the bond valence sums (BVS) 
over all Ni-S distances contained within [NiS4] tetrahedraP3 
The population is symmetrically distributed about a BVS = 
2.0, with a full-width-at-half-maximum of 0.2, as opposed to 
a non-integer value or mixed- valence distribution. 

In these analyses, there is no evidence for long-range or- 
dered magnetism. The structural models used to describe the 
high-resolution synchrotron diffraction data provide excellent 
fits to the neutron powder diffraction data [Figure [8ja,b)]. Di- 
rect subtraction of data sets collected at T = 50 K and 5 K 
[Figure [8jc)] does not reveal the appearance of any additional 
scattering, as would emerge from magnetic order. 
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FIG. 6: (a) Representative Reverse Monte Carlo simulation of 
the pair distribution function, (b) Projection of all 24000 atoms 
from the RMC super-cell (small spheres) back onto the idealized 
crystallographic unit cell (large spheres) reveals that the ThCr2Si2- 
type connectivity of KM2S2 is maintained after the simulation, as 
is the average lA/mmm crystallographic symmetry, (c) A represen- 
tative [N12S2] layer extracted from a 24000 atom RMC supercell 
(T = 300 K) illustrates the incoherent distribution of short (dark, 
^Ni-Ni < 2.52 A) and long bonds (light, ^Ni-Ni > 2.86 A). Interme- 
diate length bonds are omitted for clarity, (d) Scaled atomistic rep- 
resentations of short and long Ni-Ni bonds extracted from the T = 
300 K RMC supercell. 



B. Physical Properties 

Resistivity measurements indicate metallic behavior at all 
temperatures, as previously reported [Figure [9ja)]. 28 There 
is a discontinuity near T ~ 250 K with hysteresis, as ob- 
served from first-order phase transitions. This transition co- 
incides with a significant change in the isotropic thermal dis- 
placement parameter of the K + sublattice as inferred from 
the SXRD data [Figure |9jb)]; however, neither the S position 
[Figure [5Jc)] nor the Ni position are greatly disturbed [Fig- 
ure^d)]. 

Measurement of the linear magnetic susceptibility reveals 
only a weak temperature dependence (Figure [TO]). The mag- 
netization was measured in two ways: isothermally and at 
constant field. The constant-field magnetic susceptibility (%) 
was approximated by, % w AM /AH = [M21 — Mix] /[IT], in 
order to subtract trace ferromagnetic Ni impurities that give 
a subtle curvature to the T = 300 K magnetization, as shown 
in Figure [TTJa) (impurity concentration < 1%, undetectable 
by SXRD). The constant-field susceptibility exhibits a grad- 
ual upturn below T < 75 K. To test if this upturn results 
from a contribution of localized moments following a Bril- 
louin function, many isothermal field-dependent magnetiza- 
tion measurements were measured for T < 100 K. The curva- 
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FIG. 7: (Color online) (a) Histogram of Ni-Ni distances extracted 
from the RMC supercell. The orange (light gray) line is the result 
from a simulation of the measured data. The solid black line corre- 
sponds to a histogram of Ni-Ni distances simulated from anisotropic 
and harmonic displacements. Below is their difference, (b) The en- 
semble bond valence sums (BVS) indicate a single population of 
charged Ni species centered around a BVS = 2.0(2). The solid line is 
a fit to a single Gaussian expression. 



ture of the isothermal magnetization, M(H), pictured in Fig- 
ure [TTJa) and (b), is not well described solely by a Brillouin 
function: the magnetization for juoH > 5 T is linear, even 
down to T = 2 K. 

To extract the trace ferromagnetic impurity from the 
isothermal magnetization curves, the T = 300 K magnetiza- 
tion data were subtracted from the T < 100 K data. To account 
for the remaining curvature of the data, we performed a global 
fit of all of the magnetization data to, 



M - M300 k = ngJBj{x) + [x(T) - Xsook] H 



(1) 



where %(T) is the temperature-dependent linear susceptibility, 
%300K is the linear slope of M(H , 300 K), n is number of lo- 
calized paramagnetic spins per mol Ni, g is the gyromagnetic 
ratio, / is the total angular momentum, x = (jubH) / (ksT), H 
is the applied magnetic field, and 
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2J 



1 , f(2J+l)x 
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(2) 



From the global fit with g = 2, we extracted / = 2.0(1) and 
n = 3.2(1) x 10 -4 impurity spins per mol Ni. The values of 
%(T) are shown in Figure 10 as linear M/H. We obtain equiv- 
alent values of % by simply extracting the slope of the lin- 
ear portions of the magnetization curves [% = AM /AH; Fig- 
ure [TT|a)]. Because the values of % are very small (~ 10 -3 
emu mol Ni -1 Oe _1 ), trace impurities (Ni 2+ , / ~ 1 to 4) have 
a significant effect on the observed magnetization. 

The linear contributions to the isothermal magnetization 
follow a weak temperature dependence; fitting these val- 
ues [closed squares, Figure [TO) to the Curie-Weiss equation, 
% = C/(T — 0)+%o, allows us to calculate a lower bound to 
temperature-independent contribution of the magnetic suscep- 
tibility, %o > 4.9(5) x 10~ 4 emu mol Ni -1 Oe -1 . 

The low-temperature specific heat data (Figure [12]) reveals 
a ^-type anomaly consistent with bulk superconductivity at 
T c = 0.46(2) K. For 1.8 < T < 20 K, the total specific heat 
was modeled as, C = yT + P37" 3 + PsT 5 , to extract the elec- 
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FIG. 9: (a) Temperature dependence of electrical resistivity of 
KM2S2, measured on a poly crystalline pellet, indicates metallic 
transport and a first-order transition near T ~ 250 K. (b) The 
isotropic atomic displacement parameters for K, Ni, and S show tran- 
sitions concomitant with the K + sublattice freezing at T ~ 250 K 
(dotted line). Colored lines are guides to the eye. (c) Rietveld analy- 
sis of the SXRD data reveals that the S position remains undisturbed, 
with only a minor transition observed at T ~ 250 K. 



FIG. 8: (Color online) Neutron powder diffraction (NPD) data and 
from the NPDF instrument (90° bank) acquired at (a) T = 300 K and 
(b) T = 15 K and Rietveld analyses illustrating the absence of ad- 
ditional reflections, as would appear from long-range ordered mag- 
netism and enlargement of the lA/mmm nuclear unit cell (circles: 
data; line: fit; difference curve below), (c) Direct subtraction of NPD 
data (black) collected at T = 50 K (orange/light gray) and T = 5 K 
(blue/dark gray) on the HIPD instrument to illustrate an absence of 
additional scattering from long-range ordered magnetism. 



tronic contribution to the specific heat described by the Som- 
merfield coefficient, y = 68(1) mJ mol" 1 K~ 2 (Table |5}. The 
specific heat jump at the transition is AC e /yT c = 1.7. Small 
external magnetic fields suppress the superconducting transi- 
tion, with # c2 (0K)~0.04(l) T, obtained by fitting the observed 
field dependence of T c to a two-fluid model. The Sommerfield 
coefficient is only weakly dependent on an applied magnetic 
field [Figure [12] inset]. The small upturn in the heat-capacity 
at ju H = 14 T for T < 0.2 K is well described by Schottky 
anomalies for nuclear and impurity spins. 

Normalizing the specific heat measured to higher tempera- 
tures by T 3 [Figure [l3| reveals non-dispersive phonon con- 
tributions to the lattice heat capacity. Dispersive phonons 
should plateau when plotted as C/T 3 with decreasing 
temperature j 29 l 3Q l meanwhile the electronic contribution rises 
sharply (C e /T 3 = y/T 2 ). Therefore, the high-temperature spe- 
cific heat was fit to combination of a Debye lattice model and 
several Einstein modes describing non-disperseive, localized 
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FIG. 10: The linear magnetic susceptibility (closed squares; thin 
line as a guide to the eye), as determined by the high-field slope of 
isothermal magnetization, is nearly temperature-independent. The 
magnetic susceptibility determined by measurements at constant 
fields (open circles) reveals a gradual upturn at T ~ 75 K. 



lattice vibrations. The total heat capacity was described as: 

f x D x 4 e x dx 



C v =yT + 9to 



(e x -l) 2 



+ 



(3) 



where R is the gas constant, So is the Debye temperature, 
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FIG. 11: (a) Isothermal field-dependent magnetization curves at var- 
ious temperatures (open circles) illustrate the linear contribution at 
high magnetic field (solid lines) and increased curvature at low tem- 
peratures, (b) Global fit (thin lines) of a Brillouin function and linear 
susceptibility term (Eqn.[TJ to the M — M300 k data (open circles). 
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FIG. 12: The electronic specific heat (C p — Ci att ) reveals a bulk 
superconducting transition at T c = 0.46(2) K (thick solid line: equal 
entropy construction). Magnetic fields suppress the transition, and by 
/joH= 14 T a sharp upturn at T < 0.2 K is present and described by 
Schottky anomalies corresponding to impurity (~ 10 -6 ) and nuclear 
spins (thin solid line). Inset: applied magnetic fields have a weak 
influence on y. 
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FIG. 13: (a) The temperature dependence of the unit cell volume 
of KN12S2 shows two anomalies below T < 13 K, with an overall 
negative coefficient of thermal expansion below I < 8.9 K. (dashed 
lines guide the eye, vertical lines denote transitions), (b) The excess 
entropy below T ~ 10 K is only well described when an enhanced- 
mass yT electronic contribution is included. The Dulong-Petit limit 
(3NR/T; dotted line) is achieved near T = 300 K. (c) The specific 
heat normalized by T 3 indicates additional thermodynamic degrees 
of freedom, well described by the Einstein expression for localized 
vibrational modes, as described in the text (yT contribution: dashed 
line; Debye lattice: black dot-dashed line; Einstein modes: shaded 
peaks). 



the low-temperature curvature of the Einstein expressions and 
electronic contributions, respectively. The model provides an 
excellent fit to the experimental heat capacity data; the fit 
parameters are tabulated in Table |n| Two localized (disper- 
sionless) vibrational modes, each described by an Einstein 
expression with energies, to = 7.47(3) meV [87(4) K] and 
34(1) meV [394(12) K] are found in addition to contributions 
from Debye expression for dispersive phonon and the con- 
duction electrons. Fits with only a single Einstein mode do 
not provide description of the data, but we cannot rule out the 
possibility that there are more than two Einstein modes on the 
basis of these data. 



xd = ®d/T, s is the number of Debye oscillators, to; is the 
energy of the i th dispersionless mode, and pi is the number of 
dispersionless oscillators. By simultaneously fitting the mea- 
sured heat capacity to Eqn. [3] as C/T (in J mol -1 K -2 ) and 
as C/T 2 (in J mol -1 K -3 ), it ensures a proper weighting of 
the high-temperature curvature of the Debye expression and 



IV. DISCUSSION 

The subtle, but rich structural transitions observed in 
KM2S2 are indicative of strongly-correlated electronic 
physics. The fact that the Ni displacements disappear on cool- 
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TABLE II: Specific heat fitting parameters 



Electronic term 


Debye terms 




Einstein terms 




yCmJmol- 1 K~ 2 ) 


S D (K) s (osc.) 


hod\ (meV) 


Pi (osc.) 


h®2 (meV) 


P2 (osc.) 


67.6(1) 


219(2) 2.52(5) 


34(1) 


1.50(5) 


7.47(3) 


0.47(1) 



ing is a highly unusual observation, an indicator of an increase 
in local symmetry on cooling. To our knowledge, the only 
materials in which this has been observed are the colossal 
magnetoresistive perovskite manganitesEKl PbRuC^JSl the 
binary lead chalcogenidesp^ and the analogous compound, 
KNi 2 Se 2 P 

While many materials exhibiting a negative coefficient of 
thermal expansion (NTE) are known, in most examples, such 
as Re03,^ a NTE typically arises from the connectivity of 
rigid polyatomic units that become less flexible on cooling 
or from the increased amplitude of rigid unit modes on heat- 
ing that pulls in the structure, as with the analogy to a gui- 
tar string^. Such behavior is not expected in KM2S2 due 
to the constrained connectivity of edge-sharing [NiSzt] tetra- 
hedra within the layers. Instead, the negative thermal ex- 
pansion of KN12S2 may reflect the breaking of directional 
bonds and derealization of charge, as illustrated by caesium, 
where the addition of hydrogen results in a decrease in vol- 
ume per formula unit from Cs (111.7 A 3 /Cs) to CsH (64.8 
A 3 /CsH) ! 38 l 39 l This comparison suggests that we observe the 
formation of a more delocalized electronic state in KM2S2 at 
the lowest temperatures and analogous to the driving force be- 
hind anomalous thermal expansion in heavy fermion materials 
{e.g., URu2Si2)! 40 * 41 * These observations may also be related 
to intermediate valence compounds, YbCuAl 42 or CeA^p^ 
where a more spatially extended valence state is favored at 
lower temperatures. 

However, here, such an effect requires an involvement of 
direct Ni-Ni bonding, as the bond valence sum analysis of 
the Ni-S bonding does not show any evidence for charge dis- 
proportionation or Jahn-Teller distortions in KNi2S2.^^ In- 
stead, our results and the absence of charge order are consis- 
tent with the off-centering displacements arising from chage 
density wave (CDW) fluctuations within a manifold com- 
prised predominately of Ni d x i_ y i orbitals, as proposed for 
KNi2Se2.^ These CDW fluctuations are spatially related to 
not only the long-range ordered CDWs observed in KCu2Se2 
and SrRfi2As2 p 2 * 23 * but also the tetragonal-to-orthorhombic 
structural distortions observed in BaFe 2 As 2 46 and Fei.oiSeF^J 

Quantitative analysis of the specific heat at higher temper- 
atures reveals insight into the dynamic nature of the structural 
anomalies as charge density wave fluctuations. On cooling, 
most of the entropy of the 34(1) meV mode is released at a 
similar temperature scale to the disappearance of the CDW 
fluctuations; the entropy of the 7.47(3) meV mode is released 
at the same temperature scale as the observation of negative 
thermal expansion behavior T < 8.9 K (Figure [13]). Analysis 
of the heat capacity is consistent with both the loss of a CDW 
around T^75 K and the emergence of heavy electron physics 
at low temperature, suggesting that the two phenomena are 



interrelated. 

While a significant increase in carrier mobility concomitant 
with the depopulation of Ni displacements was observed in 
KNi2Se2,^ the transport properties of KM2S2 appear to be 
dominated by the freezing of the K + sublattice [Figure [9ja)], 
as inferred from the significant decrease in the K Ui so ob- 
tained from the SXRD data [Figure |9fb)]. Such transitions 
are usually first-order due to an associated latent heat, as per 
liquid- solid transitions pa and potassium has been shown to 
be readily mobile in the related compound, KNi2Se2 at room- 
temperature!^ This freezing process likely introduces many 
electronic scattering centers and prevents an accurate mea- 
surement of the intrinsic resistivity. 

We observe two significant thermodynamic signatures of 
many -body physics of KM2S2 in addition to the negative ther- 
mal expansion: superconductivity and an enhanced electronic 
band mass. The low-temperature specific heat data [Fig- 
ure [l2) reveals a ^-type anomaly consistent with bulk super- 
conductivity at T c = 0.46(2) K. Such a small H C 2 (compared to 
T c ) is indicative of an enhanced mass of conduction elections, 
between m# = 40m e to 80m e . This estimate of the degree 
of electronic mass enhancement depends on the assumptions 
about the size and shape of the Fermi surface. Using the ex- 
trapolated T = K upper critical field, H c2 = 0.04(1) T, the 
average zero-temperature coherence length of the supercon- 
ducting state is £(0) = [§o/(2nH c2 )} l/2 = 90(10) nmPThe 
Fermi velocity is estimated from T c and £ using the relation 
v F = (k B T£)/(0.1$h) = 3(1) x 10 4 m/s. Then, assuming 
a spherical Fermi surface, the Fermi wavevector is estimated 
from the carrier density, n = carriers per unit cell volume, us- 
ing kp = (3^7l 2 ) 1 / 3 . If there are 1.5 carriers per Ni, then kp = 
0.99(1) A -1 , while if all 33 valence electrons per formula unit 
contribute, then kp = 2.2(1) A -1 . The resulting mass enhance- 
ment is then calculated from the relation nf /m = hkp/vp, and 
ranges from nf Hc2 /m = 40 to 80. 

The specific heat in the normal state above T c also points 
to an enhanced band mass: the low-temperature Sommerfield 
coefficient, y= 68(1) mJ mol -1 K -2 , represents a significant 
electronic mass enhancement, between nf = 1 \m e to 24m e . To 
estimate the electronic band mass enhancement, we assume a 
spherical Fermi surface with either 33 valence electrons per 
formula unit, or 3 valence electrons per formula unit (Ni L5+ ). 
The carrier density, n, is given by the number of carriers, N, 
per cell volume, V, to give n = N/V, which is used in the 
calculation of the Fermi wavevector, kp = (3% 2 n) 1 ^ 3 , in order 
to estimate the densities of states at the Fermi energy, g(Ep). 
The Sommerfield coefficient from the spherical Fermi surface 
is given by y e = K 2 k^g(Ep)/3. The electronic band-mass en- 
hancement is then estimated from nf /m e = Ymeasured/Ye- 

This electronic mass enhancement, comparable to that ex- 
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FIG. 14: The temperature/magnetic-field phase diagram displays the 
complex structural and electronic states observed in KN12S2, includ- 
ing superconducting (SC), heavy mass, intermediate but non-charge 
density wave (CDW), Ni-Ni charge density wave, and high K + mo- 
bility phases. The K + sublattice melting transition is a first-order 
transition at T = 250 K with an associated two-phase region. The 
other transitions appear to be second-order, but the data are ambigu- 
ous regarding their thermodynamic nature. 



tracted from the upper critical field, represents a much larger 
mass enhancement than is generally observed in metallic 
correlated solid d 50 ! 51 ! and is comparable the mass enhance- 
ment observed in the archetypical heavy fermion compound, 
URu2Si2p^ Unlike prototypical heavy fermion materials E2HS3 
however, the Sommerfield coefficient is only weakly depen- 



dent on an applied magnetic field [Figure 12 inset], suggest- 
ing at most a minor role for magnetism or magnetic fluctu- 
ations in producing the heavy mass state. This is addition- 
ally supported by the lack of localized magnetism observed 



in magnetic susceptibility measurements (Figure 10 and Fig- 
ure [TT]) and high-flux neutron powder diffraction experiments 
(Figure [SI. 



V. CONCLUSIONS 

In short, our data show that KM2S2 exhibits a rich and 
unusual electronic and structural phase diagram below T ~ 



440 K, as summarized in Figure 14 Near room temperature, 
the K + sublattice is mobile and exhibits what appears to be a 
freezing transition near T ~ 250 K. Above T ~ 75 K, there 
are aperiodic, incoherent CDW fluctuations corresponding to 
displacements of the Ni-Ni sublattice, concomitant with com- 
plete release of the entropy from a localized vibrational mode 
at to = 34(1) meV. Below T ~ 75 K, we observe an inter- 
mediate state, which can be described as a correlated metallic 
state in which there is no CDW nor significant electronic mass 
enhancement. Upon further cooling, KM2S2 displays a sig- 
nificantly enhanced electronic mass, 11< m* /m e <24, with a 
concomitant lattice expansion [negative coefficient of thermal 
expansion] and the entropy release of a dispersion less vibra- 
tional moe of 7.47(3) meV. Below T c = 0.46(2) K there is a 
superconducting transition that is suppressed by a //^(OK) = 
0.04(1) T. Surprisingly, all of these strongly correlated behav- 
iors occur in the absence of localized magnetism. Instead, 
our results suggest that the origin of heavy electron behav- 
ior in KM2S2 lies in the hybridization o f near ly localized and 
bonded states with conduction electrons 12323 It will be inter- 
esting to establish how proximity to charge order can drive 
strongly correlated physics in this and related materials fami- 
lies. 
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